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Abstract 

The equations for various spin particles with oscillator-like interactions are discussed in 
this talk. Contents: 1. Comment on ’’The Klein-Gordon Oscillator”; 2. The Dirac oscillator 
in quaternion form; 3. The Dirac-Dowker oscillator; 4. The Weinberg oscillator; 5. Note on 
the two-body Dirac oscillator. 


1 Comment on ’’The Klein-Gordon Oscillator” 


In connection with the publications of Moshinsky et al., e. g. [1], the interest in the model with 
the j = 1/2 Hamiltonian that is linear in both momenta and coordinates has grown recently [2]. 
Analogous type of interaction has been considered for the case of j = 0 and j = 1 Duffin-Kemmer 
field [3] and for the case of j = 0 Klein-Gordon field [4]. 

In the paper [4] the operators Q, coordinate, and P, momentum, have been represented in 
n®n matrix form 

Q = m, P = vp, (Ll) 

with fj 2 = 1. The interaction in the Klein-Gordon equation has been introduced in the following 

wa Y : _ . _ 

P -» P -imjn -Q, ( 1 . 2 ) 

where for the sake of completeness f 1 is chosen by 3 <S> 3 matrix with coefficients fl t] = The 

7 matrix obeys the following anticommutation relations {7, i)} = 0, 7 2 = 1. 

The Klein-Gordon equation for '!'(<?,<), the wave function which could be expanded in two- 
component form, is then 


- = ( p 2 + m 2 q ■ 


Cl 2 ■ q + m-y trQ. + m 2 ) ^(q,t), 
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what gives the energy spectrum [4] 


E( a )N, ~ m2 = 2m (u\ N\ + U 2 N 2 + W3A3) , N U N 2 ,N 3 -0,1,2... 

E(b)Ni ~ m2 = 2m (o;i(JV| + 1) + wa(-^2 + 1) + ^ 3(^3 + 1)) • (1-4) 

However, the physical sense of implementing the matrices 17 and 7 in [4] is obscure. In this 

Section we try to attach some physical foundations to this procedure. It is well-known some 

ways to recast the Klein-Gordon equation in the Hamiltonian form e. g. [5, 6]. First of all, the 
Klein-Gordon equation could be re-written to the system of two coupled equations [6,p.98] 


dV 

dx a 


= K=., 


d=/ 

dx l 


= -/c'F, 


(1.5) 


where k = mc/fi (in the following we use the system where c = h = 1). By means of redefining 
the components they are easy to present in the matrix Hamiltonian form (cf. with [7]) 
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X2 



\X3/ 
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P 1 
P 2 


Pi 

0 

0 
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P 2 
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0 

0 \ 



( (f>\ 

0 

0 
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0 

-1 

0 
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0 

0 

0 

0 
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0 
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0 

0 ) 


Vo 

0 

0 

-l) 



Vx3/ 


(1.6) 


provided that 4> = id t V + mV, Xi — —i V> ^ Using matrices a and /?, corresponding to this 

case, and introducing interaction analogously [la] we come to the equation for upper component 


( E 2 — m 2 )0 = | p 2 + m 2 u 2 r 2 — 3mu>j (j> 


(1.7) 


what coincides with Eq. (10a) of ref. [4] in the case u>\ = u ) 2 = 0 J 3 . 

The similar formulation also originated from the Duffin-Kemmer approach. In this case the 
wave function = column( 4 >i, 4 > 2 ,Xi) is five-dimensional and its components are 4>\ = ( id t $ + 
m*)/V 2, eh = — m'I'J/v/ 2, Xi = —i Vi 1 ^ = Pi'®- It satisfies the equation 

= ( Bp + mP 0) S, B, = [A,,/?*]- (1-8) 

(our choice of 5® 5 dimension /^-matrices corresponds to ref. [3]). As shown there, the substitution 
p — f p— imuriQr leads to the equation (1.7) for both 4>\ and (j> 2 . Let us remark, in both the approach 
based on Eq. (1.6) and the Duffin-Kemmer approach, Eq. (1.8), we have the equation 

(. E 2 — m 2 )xi = ( Pi ~ imuix , )(pj 4- imux j )xj (1-9) 

1 

for down component, which seems to not be reduced to oscillator-like equation. 

Then, Sakata-Taketani approach, e. g. [5], is characterized by the equation which we write in 
the form: 

(iio) 

with Ti being the Pauli matrices. $ = eolumn(<f), y) is the two-component wave function with 
components which could be written as following: <f> = ('F + x = — ~d t ^)/ >/2. 
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From the previous experience we learned that in order to get the oscillator-like equation we need 
to do substitution with matrix which anticommutes with matrix structure of the momentum part 
of the equation. In our case the matrix which has this property is T\ matrix. Therefore, we do 
the substitution p — > p — imuT\f and come to 


E 2 ti 


p 2 + m 2 u> 2 r 2 — 3 mu + m 2 j £, 


( 1 . 11 ) 


where £ = <p + and to the analogous equation for q = <p — X — + x)- I n the process of 

calculations we convinced ourselves that the interaction Hamiltonian 

H — ^ (p 2 + m 2 (jj 2 x 2 — 3 mw) (r 3 + jY 2 ) + mr 3 (1.12) 

is the same as in [3, formula (3.9)] since rj(r 3 + *72)71 = — (r 3 + ir 2 ) and (r 3 + ?'t 2 )ti = t 3 + ?r 2 . 


2 The Dirac oscillator in quaternion form 


The quaternion (and conjugated to it) with real coefficient is defined as q = <70 4- iq\ + jq 2 + kq^, 
q = qa — iqi — jq 2 ~ kq^. The basis vectors satisfy the equations i 2 = j 2 = A: 2 = — 1 and ij = ~ji = k 
with cyclic permutations. Considering a two-component quaternionic spinor (or SL(2,H) spinor) 
one could write the free, Dirac equation as following, ref. [8c, d], 

f • d — mr^k — 0. (2-1) 

Anticommutation relations for f are given in [8d, p.222]. In Pauli representation (i — ► — \/ — Itj , 
j — ► — Ir 2 and k — ► — v^—Tt 3 ) it goes through to usual Dirac equation and its complex 
conjugate. As mentioned in [8] it is convenient to diagonalize the matrices entering in Eq. (2.1) 
using matrix 

T = 71 (1 7?)' (22) 

In such a way we come to biquaternionic formulation (g, € C): 


f dipi -I- imip ] R = 0 (2.3) 

\ d^ R + imip L = 0 , 

where ipi = ipp+, ipR = ipp~- This decomposition of ’J* into left ideals is carried out by means of the 
projection operators p± = (& 0 ±f> 3 )/2. New basis is 6 0 = l,&i = = y/—Tj, b% = \f^\k and 

b 0 = 6 0 , b Q = —b Q . Introducing interaction in the form d, — ♦ d,+r 3 Vi(x), V is the compensating field 
for this type of 5p(l, Q) transformations, and taking into account that the vectors of biquaternionic 
basis anticommute b Q bp + bpb Q = — 2 q a p, q a 0 = diag( — 1, 1, 1, 1), we come to the equations for the 
left and right spinor-quaternions in the following form: 

{E 2 - m 2 )ipi = [(p 2 + k 2 x 2 ) - 3 k - 2e ijk b k x,p j \ ip L (2.4) 

{E 2 - m 2 )$ R = [(p 2 -I- k 2 x 2 ) + Zk + 2 e ijk b k x i p j ] \p' R (2.5) 

if we choose Vi(x) ■= kx { . Eqs. (2.4) and (2.5) are the Dirac oscillator equations in the Pauli 

rep, b k — ► T k . Analogous equations for xpR and could be obtained from (2.3) if one choose the 
opposite signs at the mass terms. 
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3 The Dirac-Dowker oscillator 

In this Section we start from the equation for any spin given by Dirac [9] in the form written down 
by Corson, ref. [6,p.l54], (here we use Corson’s notation) 

f d AB v B (k+ + \,l-\)-m (y^Y /2 v A (l)ip(k,l) = 0 

{ ( 3 - 1 ) 

l d Ag v A (l)xp(k,l) + m (iHti) 1 ' 2 v B (k + \)^{k + J,/ - \) = °, 

where v A and are the rectangular spinor-matrices of 2k rows and 2k + 1 columns (see, e. g., 
section 17b of ref. [6] for the details). The wave function ip(k,l) belongs to the ( k,l ) representation 
of the homogeneous Lorentz group. The choice 1 = 1/2 and k = j — 1/2, j is the spin of a particle, 
permits one to reduce a! number of subsidiary conditions. Moreover, the equations (3.1) are shown 
by Dowker [10] to recast to the matrix form which is similar to the well-known Dirac equation for 
j = 1/2 particle 


' = mT 

I (3.2) 

a^d^T = —md>. 


The 4 j- component function $ could be identified with the wave function in ( j , 0) © (j — 1 , 0) 
representation. Then, T, which also has 4 j components, is written down 


T = 



v aU ~ 

u A (j) ® 1^(5) ) 




(3.3) 


and it belongs to ( j — 1/2, 1/2) representation. The matrices a* 4 and = a h obeys all the 
algebraic relations of the Pauli matrices = g M ", except for completeness, i 

Defining p ^ = —id^ and the analogs of 7- matrices as following: 


Y = 


(° 

[ia * 


— ia * \ 

0 ) 


(3.4) 


the set of equations (3.2) is written down to the form of the Dirac equation 

(P»Y ~ m) = °. (3.5) 

However, let us not forget that $ and T are 2-spinors only in the case of j = 1/2. 

In the case of spin j = 1/2 it is well-known the set of 7- matrices is defined up to the unitary 
transformation and Eq. (3.5) could be recast to the Hamiltonian form given by Dirac (with a* 
and P matrices) by means of the unitary matrix. It is easy to carry out the same procedure 
(a k = <S7°7*\S -1 and /? = S7°S~ 1 ) for 7- matrices, Eq. (3.4), and functions of arbitrary spin 
(\t = <S _1 $). For our aims it is convenient to chose the unitary matrix as following: 

g _ 1 i^4j<s>4j \ ^2 

H 4j '. 'Aj A 4j®4j / 


336 



After standard substitution p — ► p — imu>i°r we obtain 

E<p = —i [ao(ap) + ima;(af)] v + (3-7) 

Ev = i [a 0 (ap) — imu^ar)] — ma 0 f. (3.8) 

Since it follows from thy anticommutation relations that a^ao = ao a i we have the equations which 
coincide with Eq. (8) of ref. [la] or Eqs. (3.6) and (3.12) of ref. [Id] except for — * a h , i. e. 

their explicit forms, 

(E 2 — m 2 )<t> = [p 2 + muf 2 - 3a 0 rau> — mwaoa^a^V'Vj] 4> (3-9) 

(E 2 + m 2 ) v = [p 2 + mujr 2 + 3a<)ma> + muoiQO^dj^r' Vj] V- (3.10) 

Thus, we convinced ourselves that we got the same oscillator-like interaction and the similar 
spectrum as for the case of j = 1/2 particles in [la]. 


4 The Weinberg oscillator 


The principal equation of 2(2 j + 1)- component 


approach [11] in the case of spin j = 1 is 


(-WWV + M 2 )¥i='\x) = 0, 


(4.1) 


with j a 0 being 6 ® 6 covariantly defined matrices. The j — l Hamiltonian has been given in refs. 
[llb,c]: 


H = 


‘IE 2 

2 E 2 - M 2 


(ap) + 0 


E - 


2E 


(«P) 2 1 i 


2 E 2 - A/ 2 


(4.2) 


where 



0 J 


(Si are the spin matrices for a vector particle). 

In general, the upper and down components of 6- component wave function do not uncouple 
neither under the interaction p — ► p — imu)(3r nor under 75,/ii/ u /» T V However, if we introduce 


the Dirac oscillator interaction so that the conditions of the longitudity of $ = column^, \i) 
respective to r, i. e. rx</> = 0, rx\ = 0 are fulfilled, we come to the equations more simple 

(2 E 2 - M 2 )£ = E(Sp)ri + [(Sp) - fc(5i 1 ')] (Sp)£, (4.3) 

E(SM = [(Sp) + k(Sf)] (Sp)r] (4.4) 

(f = (j) - x, rj = <p + x), which could be uncoupled to the following form ( k = imu ) 

(Sp)(E 2 - M 2 )(Sp)i = ( Sp) [p 2 + mu 2 r 2 + 3 mu + AmuS[r x p]] (Sp)£ (4.5) 

(Sp)(E 2 - M 2 )(Sp)r] = ( Sp) [p 2 + mu) 2 f 2 - 3 mu - 4muS[r x p]] (5p)p - 

—im<jj(2E 2 - M 2 )(Sr)(Sp)ii (4.6) 


These equations can be considered as the extension of the equations with Dirac oscillator 
interaction to the j = 1 case, for the components (Sp)£ and ( Sp)r ]. However, remark that one has 
the additional spin-orbit term acting as earlier at rj. 
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5 Note on the two-body Dirac oscillator 

The two-body Dirac Hamiltonian with oscillator-like interaction is given by (see, e. g., ref. [lc]) 


Hip == 


|^(Si + a 2 )-P + ^(<Ji ~ A*) • P~ “ "2) • rB + m(pi + AO 


ip. (5.1) 


In the c.m.s. it is possible to equate P = 0. The matrices are given by the direct products 


- ( 0 
ai = U 

A 

0 

MY 

0 ) 

$■202/ 

, «2 = 

■( 

^2®2 
0 1 

iL) ® 1 

rs 

K°2 

#2 \ 
or 

(5.2) 


B = 

A ® As = | 

( H2®2 
^ 0 

0 1 
— fl2®2^ 

I <s> 

/ l2®2 
V 0 

0 1 

~ H2®2/ 

1 


(5.3) 


r 5 

= 7 ? ® 72 = 

-{ 0 
\l2®2 

ti-202 > 
0 > 

) <8> 

( 0 
\ l2®2 

H2®2'\ 

0 ) ' 



(5.4) 


Now we apply the same procedure like that was used for transformation the Bargmann-Wigner 
equation to the Proca equations. The 16- component wave function of the two-body Dirac equation 
could be expanded on the complete set of matrices: ( 7 ^C), (a^C) and C, ( 7 5 C), and (y 5 y fi C). 
We consider the system multiplied by C, the matrix of charge conjugation, in order to trace for 
the symmetric properties under oscillator-like potentials. The wave function is decomposed in 
symmetric and antisymmetric parts using the above-mentioned complete system of matrices: 


V’{q/ 3} — 7Qf)Crj/3 J 4/j T v 

*P{a{f\ = C Q p<p + 1 0^7)13$ A 

In such a way we obtain the set of equations: 

EAq — 0, EAo — — 2m0, E<p = 2iV2(pi — ir')F'° 

E4> = -2mAo + \f2.€ijk{pi + 

EA‘ = — *V 2€ ijk (pj T ir j )A k 
EA' = 4 imF 0i + iV2c ijk (pj ± if j )A k 
EF 0i = -2 imA' -I- i\f2{pi + ir')<p 

EFj k = -j=€ ijk (pi - ir')4> 


(5.5) 

(5.6) 

(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 


Let us mention that for another type of Dirac oscillator-like interaction ~ (c?i — a 2 )BT^ the only 
changes are the sign changes at the term if in Eqs. (5.9) and (5.10) of the above system. The two- 
body Dirac oscillator equations in the form (5.7)-(5.12) could be uncoupled on the set containing 
only functions <p , 4> and A^ and the another one containing only A^ and F^„: 


1 ) 


( E 2 — 8 m 2 )<p = 4 (pi — if')(pi + ir')<p 


_1_ (7l Gmeijkftpj 
E\{ 0 


)) 


(5.13) 



(5.15) 


EAq = —2 m4> 

(E 2 - 8m 2 ) A' = 2 (pj T ir j )(pi ± if i )A j - 2 (pj T ir j )(Pj ± ir’M* + 


+|{( 16me f'' A )}« (5.16) 

2) EA 0 = 0 (5.17) 

(. E 2 - 8 m 2 )F 0 ' = 4 {pi + iP){pj - ir j )F°> - 

—4 i~p;(Pj ± ir*)(pi ir* )A* + 4«— (p, ± i^){Pj 4- ir^)A l (5.18) 

E E 

E 2 A‘ = 2(;7j ± ir J )(/7, ^ ,7").4 J - 2 (pj ± ir’)(pj =F ir>)A‘ + 4imEF 0i 

(£ 2 - 4m 2 )fl‘ = - iF)(p, + ir‘)F m ". (5.19) 

I 

1 

This fact proves the Dirac oscillator interaction, like the case of introduction of electrodynamic 
interaction in the Proca or the Bargmann-Wigner equations, does not mix 5=1 and 5 = 0 states. 
Next, the interaction term of the following form: 


ynl 


1 dV{r)/dr 
r 1 — [V^r)] 2 


(<*! - a 2 )BT 5 f 


(5.20) 


has been deduced [12] from the equation of Relativistic Quantum Constraint Dynamics (RQCR) 
or N- particle Barut equation. In [12] it proved to lead to the Dirac oscillator-like interactions 
provided that the definite choice of the function V^r). In connection with that let us remark the 
curious behavior of the another potential V(r) which has been proposed in ref. [13b, c]: 


V(r) 


2 coth (rmn) 
47IT 



coth («r) 
4xr 


(5.21) 


It could be deduced from the one-boson exchange quasipotential V(p,k,E ) = —g 2 {p — k) 2 by 
means of the transformation into the relativistic configurational representation (RCR) using the 

complete set of Shapiro plane-wave functions: £(A,f) = (Ao — A n/m) _1_,rm , Ao = \/a 2 + m 2 , 
n = f/|f|. 

In the case of the quasipotential (5.21) the interaction term V, Eq. (5.20), has the different 
asymptotic behavior in three regions (<? 2 /(47r) = 1). Namely, 


V ,n< ~ 7 1 (<?i - o 2 )Br 5 r ~ 

r(r z — 1) 

' (l/r 3 )(c*i — a 2 )Br 5 r, if r >> £ and r>l 

~ < 

k -(l/r)(o! - a 2 )Br 5 r, if J « r < 1, 
in the infrared region (r >> large distances); and 


(5.22) 


V int „ _ 2k (c?, _ a 2 )Br 5 r, if r « -, (5.23) 

in the ultraviolet region (small distances). In one of the regions one has the Dirac oscillator- like 
behavior. 
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